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1.- Introduction. 



One of the most transparent pictures of the Scrodinger's eigenvalue equation in 
1-dimension is obtained by reinterpreting the space coordinate x as the "time" 
{x — t)^ and the values of the wave function ^{x) and its derivative i)'{x) as the 
position and momentum of a classical point particle [|[ ^, U, |[ |^ . The resulting 
classical trajectory obeys Hamilton equations which can be further simplified 
by introducing the Priifer's angle on the Poincare plane The integration 
of the 1-st order "angular equation" then turns out one of the most efficient 
methods to determine the energy spectra [0, D . In this note we show that the 
'classical algorithm' works even better for the Helmholtz problem ^, |7[ ||; in 
addition, it exhibits a new aspect of the Helmholtz spectra. 



2.- The Helmholtz spectrum. 



Given an interval [a, 6] C R and a real function (j) : [a, b] — > R, the Helmholtz 
equation is a particular case of a Sturm-Liouville problem (see 0, [sj): 

t^ + \4>{x)i;{x)^0 (2.1) 

and the Helmholtz spectrum is defined as the set of all A G R for which (2.1) 
has a non trivial solution il) vanishing on both ends: 

V'(a) = = (2.2) 

The 'trajectory picture' ||l|, |^ applies quite naturally to the Helmholtz eigen- 
problem, and the spectral condition is much simpler than in the Schrodinger's 
case. Let ipi^) be any real solution of (2.1). Define: x — t, ip — q, ip' — p; eq. 
(2.1) then becomes: 



dq 

-dt = ^ 

% - -m)q (2.3) 



with the 'spectral condition': 



q{a) = q{h) ^ (2.4) 

Observe, that (2.3) are the canonical equations for a classical point-particle 
with the Hamiltonian: 

iJ(t) = ^ + A0(t)^ (2.5) 
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Each integral trajectory of (2.3) 'portraits' a wave function ip{x) and its first 
derivative. We shall further assume that ^ (j){t) > is a piece-wise continuous 
function in [a, b]. The attractive oscillator force in (2.3) then drives the classical 
trajectory clockwise around the origin of the phase plane (Fig. 1) . When A grows 
to +00, the end point {q{b),p{b)) must cross the p— axis an infinite number of 
times (Fig. 2). Each time this happens, A belongs to the point spectrum of 
(2.1). In this way the Helmholtz eigenproblem (2.1) has a discrete spectrum in 
the form of an infinite ladder of positive eigenvalues determined by the classical 
trajectories (compare the theorems |p|, |lO|). 

The calculation of the spectral values can be carried out by introducing the 
polar variables on the phase plane: p — pcosa, q — psina. The variable a then 
separates and (2.3) reads 



a = cos^ a + \<f>{t) sin^ a (2-6) 
P 1 - 



P 2 
with the initial and spectral conditions: 



sin 2a (2.7) 



initial condition ck(0) — (2-8) 

spectral condition Aa = a(b) — m: (n — 0, ±1, ±2, . . .) (2.9) 

Some exact solutions of (2.1-2) were recently studied by Goyal et al to check 
the validity of the WKB and Langer approximations Q. For curiosity, we 
applied the angular method to the same collection of cases. Eq. (2.6) was inte- 
grated numerically for a net of A-values by applying the Runge-Kutta method 
and the condition (2.9) selected the eigenvalues A quoted in Table 1. 





[a,b] 


Eigenvalues A 


Amplification 






'Exact' 


Approximate 


Angular 


constant a 


[t + ny 


[0,n] 


0.00174401 


0.00188559 (WKB) 


0.00174401 


-0.543046 






0.00734865 


0.00754235 (WKB) 


0.00734843 


0.51828 






0.0167524 


0.0169703 (WKB) 


0.0167517 


-0.510087 


l.le* - 1 


[0,1] 


11.145735 


11.145736 (LAN) 


11.145735 


-0.591942 






47.301484 


47.301493 (LAN) 


43.301493 


0.54752 






108.428802 


108.428616 (LAN) 


108.428629 


-0.52727 


(1 + sin27rt)^ 


[0,1] 






5.3347146 

34.1068933 

86.8947093 


-2.35707 

3.1455 

-3.42553 



Table 1: Helmholtz eigenvalues known in the literature compared with the 
results of the angular method. In addition we have examined (f)(t) = A(l -|- 
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sin27rt)^ in [0,1]. The last column gives the 'optical information', disregarded 
in the literature. 

As one can see, for the low eigenvalues the ranking is won by the almost ig- 
nored angular algorithm (both, in number of digits and simplicity of operation). 
Its additional advantage is that it uncovers the links of the Hclmholtz spectrum 
with the techniques of generating the squeezed states. 



3.- The trajectory optics. 



In a typical Sturm-Liouville problem, one looks for a solution of a differential 
equation which fulfills a definite boundary condition . The second independent 
solution is dimissed as unphysical. However, we shall show that for (2.1) the 
second solution too carries an important message. As eqs. (2.3) are linear, the 
corresponding q(t), p{t) depend linearly on the initial values q{a), p{a): 

'^W^=.(t,a)^K) (3.1) 



p{t)) "^'''"^Vp(«) 

where u{t, a) is the simplectic evolution matrix |ll| , ^ (also: the transfer 
matrix). The existence of a non-trivial solution of (2.2) with q{a) = q{b) = 
means that the evolution matrix u(a, b) in the interval [a, h] maps: 



and since w(a, h) is simplectic, it must be of the form: 



(3.2) 



"(&,a)=(^ a,77eR,a^O (3.3) 

Once this is known, one can deduce not only the spectral orbit of (2.3), 
but also the form of other trajectories. Consider therefore the congruence of 
trajectories which diverge at t = a from the initial point {qa^p), where q^ is 
fixed and p varies. If A is an eigenvalue of (2.1) and u(6, a) is of the form (3.3), 
then &i t = h: 



q{b)\ _(a \ (qo\ _ ( 



p{h) \r, 1/a \p \p' 



<yqo 



(3.4) 



i.e., all orbits departing from the same qo (with different momenta p) meet 
at aqo (with momenta p' = aqo + p/a). Thus, the eigenvalues A of (2.1) are 
the real numbers for which the time evolution (2.3) focuses any congruence of 
classical orbits diverging from any common initial qo- The optical constants ct, 
r] can be read from the second (disregarded) solution of (2.1). 
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Note that the focusing, in general, is not the scale transformation p2[ , (as 
the matrix (3.3) is triangular), though it is an example of squeezing ^4|, |l^ . 
The scale transformation (77 = 0) is an exceptional effect which can happen only 
for some (/)'s (see p^), whereas the focusing is a typical phenomenon which must 
occur for any positive (j){x), on an infinite ladder of A's. The converse problem is 
as well interesting. Instead of fixing an interval and looking for A, fix any A > 
and P(t) > and consider an expanding time interval [0,t]. If (f)(t) mantains 
some positive average in [0, 00), the final point (q(T),p(r)), as r — > cx3 must cross 
infinity of times the p -axis, each time generating our focusing effect no matter 
the exact form of P(t). Should P(t) be a random variable performing errating 
oscillations, the motions (2.3) will be chaotic too; yet, for any congruence of 
trajectories, a clean 'echo effect' of the initial configuration (an amplified or 
reduced image) must emerge now and then against the chaotic forces. Note also 
that the effect must have a quantum equivalent. 



4.- The Quantum Images. 



Consider the Hamiltonian H{t) in form (2.5), where q and p are position 
and momentum operators. Since (2.5) is cuadratic, the quantum evolution 
follows faithfully the classical motions. The techniques of the classical "piloting 
trajectory" gives specially good results for the Gaussian wave packets 0. 
We shall show, that it provides also an explicit expression for the evolution 
operator U{t, a): 

^^^^^ = -iH{t)U[t, a); U{a, a) = 1 (4.1) 
With this end define: 

W{t) ^ e"'^'^""U{t,a) (4.2) 

where Hq = ^{p'^ + q^) is the traditional oscillator Hamiltonian and a{t) 
is a real, differentiable function. Due to the well known transformation rules: 
^taHo q^-iaHo _ g COS tt -|-psina, Q^o'Ho pg-taHo p^osa — gsina, the operator 
W{t) fulfills: 



dW ... 2 ^ , / \ • 9 nP^ 

— I [(a — COS a— A0(i)sm a)— + 

+ (a-sin^a- A(?!)(i)cos^a)y + (4.3) 
+ (1 - X(j){t)) sin a cos a *^^ ~l ]W (t) 
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Notice, that the time dependent term coupled to corresponds exactly 
to the Priifer equation. Thus, if a{t) is the Priifer angle (2.6), then (4.3) reduces 
to: 



dW . „ qp + pq , 

-dt- ^ ^[^^Bm2a^ + 

q^ 

+ {1- X(l){t)) COS 2a^]W{t) (4.4) 
and integrates explicitcly in terms of the Priifer's radious p: 

T4^(t)=e''"(f^)^e^''(*)4 (4.5) 

Inserting (4.5) into (4.4) and remembering that the first exponent in (4.4) 
is the scale operation, one gets: 

V{t)= j\\-\4>mP^fcos2adt (4.6) 

and returning to (4.2) one has an explicit formula for the unitary operator 
U{t,a) with three exponential terms piloted by Priifer variables (2.6-7): 

Uit,a) = e-»W^V^'"(^)^e'''(*)^ (4.7) 

If in addition the number A is the eigenvalue of the Sturm-Liouville problem 
(2.1) ^ a{b) — nn {n = 1,2,...), then the first factor in (4.7) becomes P" 
(where P = parity) and 

{/(6,a) = P"e'>"(fS)^e'''^ (4.8) 
transforms any initial wave packet ijjix) — 'ijj{x, a) into: 

^u{x) = [C/(6,a)V](x) = -^e'^'f^iPix/a) (4.9) 

where a — ±p{a)/ p{b). The probability density p = \ip{x)\'^ is simultaneously 
mapped into its own diffractionless image: 

pu = = T^^p{x/<J) (4.10) 

It thus turns out that the authors looking for the squeezing parameters were 
indeed solving the XlX-century spectral problem ! The relative simplicity of the 
three term formula (4.7) (compare |l|, |l|, ^ |l| ) is due to the " bottom struc- 
ture" (2.6-9). Its additional implication concerns the " potential kicks" 6{t)V{q) 
considered recently as a method of controlling and cooling quantum systems 
While the "kicks" are mathematically simple, they are overidealized ( no 
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(5-pulses available in laboratories!). Notice, however, that if A is an eigenvalue 
of (2.1) and if 4>{t) is symmetric (f>{—t) = (j){t) in a symmetric operation interval 
[—a, a], then the eigenvector ip (and consistenly the phase trajectory q{t)) is 
even/odd in [—a, a], implying the trivial squeezing a = ±1, and the only non 

trivial part of (4.7-8) is P"e"'"2". For n = 21, this is the unitary operator de- 
scribing the results of a sharp kick of the oscillator potential ~riS{t)q'^ /2. Thus, 
the spectrum of (2.1) tells in addition how to "imitate softly" the effects of S- 
kicks converting the results on kicked oscillators ijl^ into a branch of realistic 
technology. 



5.- Magnetic Images. 



The results (3.4), (4.3) have a simple laboratory model. Consider a non- 
relativistic charged particle in 3-dim, in a homogeneous time-dependent mag- 
netic field B(t) = B{t)n, with a vector potential: 

A(x,t) = -ixxB(i) (5.1) 

where n is a constant unit vector and B(t) = Bj{t/T) is the magnetic field 
intensity, B defining 'the amplitude' and 7 the 'shape' of the magnetic pulses. 
If n point the z-axis, the Schrodinger's wave packet obeys the Hamiltonian: 



H{t) = ^pl-f3j{T)Ms+ (5.2) 
which we have written in the dimensionless coordinates: r = t/T, {qi, q2, qs) 



^^m/Th{x,y,z), (pi,P2,P3) = i\/mh/T) ^{pj:,py,pz); all physical information 
now contained in the 'shape function' 7(t) and in the 'dimensionless amplitude': 

P='^ (5.3) 
Zmc 

In (5.2) the term means the free propagation along the z-axis, the 
— /97(r)M3 generates rotations in the qi , (72-plane. The last term (a 2-dimensional 
oscillator), splits into two twin 1-dimensional oscillators of the form: 

/i(r)=pV2 + /3V(^)'zV2 (5.4) 

where {q,p) means either {qi,pi) or ((72,^2). Quite evidently, (5.4) is a 
model of (2.5) with A = ^^^^ ^ ^2^^-^ > ^(-^^ ^ ^^le spectrum 
of the Helmholtz equation (2.1) forms a ladder of eigenvalues, and if A = 
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is one of them, both oscillators yield the simultaneous focus of the trajectory 
congruencies (in classical case) or probability patterns (in quantum case). To 
check this, we have asked the computer to draw congruencies of classical tra- 
jectories diverging from a common initial point, with a common initial Vz, for 
\ = 0^ coinciding with one of the eigenvalues of (2.1) and with 0(r) = 0^{t) 
taken from Table 1. In all cases we have obtained quite precise optical images 
on planes perpendicular to the solenoid axis (see Fig. 3). The dependence of 
the phenomenon on the Sturm-Liouville spectrum is illustrated by sequences of 
amplified focuses corresponding exactly to the nodal points of the wave func- 
tion (2.1-2). Henceforth, the integration of the first order differential equation 
of Priifer (for one variable) offers the simplest algorithm to predict the optical 
echo (squeezing) for charged particle beams. 

The " reduced images" might imply new schemes to produce narrow particle 
beams. The "amplficd images" are no less interesting. They can certainly be 
used to separate the wave packet components in the quantum state engineering 
[|8[ 1^, HI). The magnetic pulses B{t) with ^(t) = -('^{t) symmetric in the 
operation interval and with A = in the spectrum of (2.1), provide a magnetic 
model for systems under the influence of oscillator kicks. 

Acknowledgements. The support of CONACYT, Mexico, is acknowled- 
ged. 
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Figure Captions 



Fig. 1.- The orbit of the classical oscilator (2.5) for <p{t) = sin^ t, \ = 10. 
The points of intersection with the p-axis determine the time moments ti = 
1.43610, t2 = 2.53311, ^3 = 4.46714, i4 = 5.51549, and the amplifica- 
tions cti = -0.52322, 0-2 = 0.65528, 0-3 = -0.23823, 0-4 = 0.27464, ... of 
the optical images. The similar phenomenon would be generated by a chaotic 



Fig. 2.- The orbits of (2.3) and the Helmholtz spectrum of (2.1). If the 

end point of the orbit with q(a)=0 places itself on the p-axis, A is an eigenvalue 
of the Helmholtz equation (1). Above: a sequence of orbits reaching q(b) = 
after multiple intersections with the p-axis illustrates the sequence of Helmholtz 
eigenstates for (f){t) ~ (1 +sin27rt)^, [a,h] = [0,1]. The corresponding optical 
constants are reported in Table 1. 



Fig. 3.- The optical effect of the Helmholtz spectrum observed in the 
solenoid for sinusoidally oscilating 7(t) = 1 + sin27rT. a) When A = = 
5.3347146 is an eigenvalue of (2.1), the congruence of particle trajectories forms 
an exact optical focus at t = 1 {t = T) (in quantum case, an amplified diffrac- 
tionless image) with a = —2.35707. The integration of the first order differential 
equation of Priifer for one angular variable turns out the simplest method to 
predict the effect, b) For X = (3'^ = 86.8947093 the congruence forms 3 increas- 
ingly amplified focuses (images of the initial point) in exact correspondence to 
the 3 nodal points of the wave function tp{x). 
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